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Abstract 

Using the operator method, the Green's functions of the Dirac and Klein-Gordon equations in 
the Coulomb potential —Za/r are derived for the arbitrary space dimensionality d. Nonrelativistic 
and quasiclassical asymptotics of these Green's functions are considered in detail. 
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I. INTRODUCTION 



When calculating the amplitudes and probabilities of QED processes in the field of heavy 
atoms, one should be aware that the parameter Za iZ is the atomic charge number and a 
is the fine-structure constant) is not small in this case. The effect of higher orders in Za 
can change the Born result by several times. Therefore, it is often required to calculate the 
probabilities of QED processes in such a strong field exactly in Za. The most convenient way 
to perform this calculation is the use of the exact Green's functions of the Dirac equation (or 
Klein-Gordon equation) for a charged particle in the field (Furry representation). Therefore, 
derivation of the Green's functions for specific field configurations is very important for 
applications. For the case of the Coulomb potential, a convenient integral representation 
of the Green's function G(r,r'|£:) was derived in Ref. [l| using the 0(2,1) algebra. The 
representation obtained is valid in the whole complex plane of the energy e and does not 
contain contour integrals. Another integral representation for the Green's function in a 
Coulomb field was derived in Ref. {2! by using the explicit form of the expansion of (?(?", r'\e) 
with respect to the eigenf unctions of the corresponding wave equation. The representation 
of the Green's function obtained in Ref. {2! contains a contour integral, which complicates 
its use in applications. 

In the calculation of the loop diagrams, it is often required to regularize the divergent 
integrals. One of the most convenient methods of the regularization is the dimensional 
regularization. In order to use the dimensional regularization within the approach based 
on the Furry representation, it is necessary to derive the exact Green's function in the 
Coulomb field in arbitrary, not necessarily integer, space dimensionality d (the space-time 
dimensionality is (i + 1). In the present paper, we solve this problem by generalizing the 
Green's function, obtained in Ref. [l| for c? = 3, to arbitrary d. Our derivation closely 
follows the path of derivation in Ref. In contrast to the conventional approach, the 
operator method used in Ref. {l| and in the present paper does not require the knowledge 
of the explicit form of the wave functions, which is difficult to define for non-integer d. In 
order to fix unambiguously the explicit form of the Green's function for arbitrary (i, we use 
only the commutative and anticommutative relations for the operators and 7-matrices . 
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II. CALCULATION OF THE GREEN'S FUNCTION 

Following Ref. we represent the Green's function in the Coulomb potential U{r) = 
—Za/r (the system of units ^ = c = 1 is used), 

G{r,r'\e) = -. S{r-r'), V = + Za/r) - jp , (1) 

V — m + iO 

as follows 



G{r,r'\e) = {V + m)D{r,r'\e) 
D{r,r'\e) = —i ds exp <2iZa£s — is 



' 2 2 K 

rp„ + K r H 

r 



5{r-r') 

5{n-n) 



K = \/m? — e"^ , Pr = — , , -^r^'^ ^^^"^ , n = r/r, n' = r'/r' 

K = f -iZacxn- {Zaf + ^{d-l){d~?,) , a = -f^-f . (2) 

Here —l^ is the angular part of Laplacian determined by 

^ = ;:^9y-%.-^f. (3) 

and 7-matrices obey the usual relation ^'^7'' + ^'^^m = gi^i^_ 
Then we represent the angular part of the 5-function as 

6in-n') = Y,Pxiri,n'), (4) 

A 

where the projection operators Px{n,n') satisfy the relations 

KPx{n,n') = X{X + l)Px{n,n'), 

J dn'Pxin, n')Py{n', n") = S^yPxin, n") . (5) 

Since the operator K contains only one matrix operator an, the matrix structure of the 
projection operator Px{n,n') is given by the linear combination of the unit matrix / and 
matrices an, an', and (an) (an'). All other matrices, such as (an) {an') {an), can be 
reduced to the four above mentioned matrices using the anticommutative relations. Taking 
this property into account, we search the projection operators Px{n,n') in the form 

Px{n,n') = aiA+(n)A+(n') + a2A+(n)A_(n') 

+a3A_(n)A+(n') + a4A_(n)A_(n') , 

A±(n) = ^(l±an), (6) 



where Oj are some functions of x = nn' . From Eqs. (|5]), we obtain 

oi = /3(A + iZa)Bn{x) , a2 = = (3 {n + u + 1/2) An{x) , 04 = /3(A — iZa)Bn{x) , 
A = ±7, 7= V(^ + ^ + l/2f -(W, /5 = ^^' 

-^1, (7) 

where C^{x) is the Gegenbauer polynomial, and n = 0,1,2,... is integer number. This 
integer number appears from the requirement that the functions Oj have no singularities at 
X = 1. The result ([7]) for was obtained with the help of the identity 

/(I + „n' + nn" . n'n")B„(nn')B„(n'n") dn' ^ ..(1 + n„")B„(nn") . 

[ , 271^^/2 271*^+1 

Vl^= dn= ——- = — , (8) 

J T{d/2) r(z/ + l)' 

Finally we obtain for projection operator 
Px{n,n') = ^ 



\[1 + {(xn){a.n')] + iZa{(y.n + an') Bn{x) 
+ {n + u + 1/2) [1 - (an)(an')]v4„(x)| . (9) 



For d = 3, this projection operator coincides with that found in Ref. j 

Note that the functions An{x) and Bn{x) have non-singular limit at z/ — > (or (i — )■ 2) 
^ sin((n + l)0) + sin(n0) ^.^ _ sin((n + 1)0) - sin(n0) 



I'-i'O sin(p sin0 

with = arccosx. 

In order to complete the calculation of the function D(r,r'\e), Eq. (|2]), it is necessary 
to find the result of the action of the operator exp{— is [rpl + H?r + A(A + 1)/t]} on the 
function 6{r — r')/r^'^. This can be done exactly in the same way as in Ref. jlj]. The method 
of Ref. jli is based on the commutator relations of the operators Ti = \[rp1 + A(A + l)/r], 
T2 = rpr, and T3 = r which coincide with those of the 0(2, 1) algebra generators (some other 
examples of applying the 0(2,1) algebra in a Coulomb field can be found in Refs.0, sl). 
The only difference between the case of arbitrary d and d = 3 is the value of the parameter 
5 in the equation Tir^ = 0. For arbitrary d, we have 

(5 = A + ^— atA>0, (5=|A| + ^— atA<0. (10) 



The final result for the function D{r^r'\e) in Eq.([2]) reads 

iViv + l) 



D(r,r'\e) 



n=0 



/ ds exp[2iZ«£s + m(r + r') cot(Ks) — 2717] 



X 



y 



J'2^{y)[l + (an)(an')]i?n(a;) + iZaJ2'y{y){oLn + Q;n')i?„,(a;) 



+ (n + z/ + 1/2) J27(l/)[1 - (an)(an')]A„(a;) j> , y = ^^7:^ , 



(11) 

smi^KSj 

where ^27(1/) is the Bessel function, An{x), Bn{x), v and 7 are defined in Eq.(I7j). The 
corresponding result for the Coulomb Green's function of the Dirac equation in d space 
dimension has the form 



G{r,r'\e) 



T = [1 + {an){an')]^^j!^^{y){-f^e + m) - iZaJ2^{y)-f'^Kcot{Ks)^Bn{x) 

[1 — {ocn){an')]{'y^e + m) — k cot(/ts)(7n — 772') J2-yiy) (n + z/ + 1/2) v4„(x) 

{-fn + -fn')J2,{y)B4x). (12) 

For d = 3 this result coincides with the corresponding result of Ref.[l|. The function 
G{r, r'\ e) has, in the complex plane e, cuts along the real axis from —00 to —m and from 
m to 00, which correspond to the continuous spectrum, and has also simple poles in the 
interval (0,m) for an attractive field and in the interval (— m, 0) for a repulsive field. The 
integral representation (|T2l) is valid for any e which belongs to the domain Ree < 0, Ime: < 
or Ree > 0, Ime > 0. If Ree < 0, Ime: > or Ree > 0, Ime < 0, then it is necessary to 
perform the integration over s in Ea. ffT2l) from zero to —00. 

n 

For real e in the interval — m < e < m we obtain (cf. Ref. [l[) 

7r/2 



00 „ 

/ dsexp[2iZaes + iK{r + r') cot{Ks) — i7i'j]T 

n=0 { 



+ 



tK^ir - r ) . ^ Q 
o • 2/ N + ^mZa-f 
2sm [Ks) 



G{r,r'\e) 



r Z/ + 1 



Ak sin[7r(Zae//t — 7)]?? 
X exp [—2iZaes //t + m(r + r') tan s\ T 

T = \l + {ct.n){oLn')] —J2 {v){'j^6 + m)—iZaJ2'y{v)'y^Ktans Bn{x 



+ 

+ 



[1 — (Q;n)(an')](7°e + m) — k tan 5(7/1 — 7n') J2'y{v) {n + u + 1/2) An{x) 

2K\/rr' 



iti^ir — r' 
2 cos^ s 



+ imZa'y 



(7n + 7n') J27(^^)5n(a;) 



coss 



(13) 



The denominator in Eq. ( !T3|) is zero at points Zaej k — 7 = for any integer k. However, 
the integral over s also vanishes for negative k at these points (see Ref. so that the 
expression (IT^ has poles only for A; = 0, 1, 2 . . .. Taking into account that 7 is positive, we 
find that the simple poles corresponding to the discrete spectrum are at the points 

m signZ 



i + (^r 



(14) 



The maximal value of Z when all the results obtained are applicable is determined by the 
relation (Za)max = [d — l)/2 (see the definition of 7 in Eq.(I7j)). 

For completeness, we also present the final result for the Coulomb Green's function of 
the Klein-Gordon equation, 

n ( '\ \ r(z/ + l) ^n + z/ 

^ ' n=0 

00 

— - — - exp[2iZ«es + m(r + r') cot(/ts) — i'iTfj]J2u(y) , 
sin(Ks) 



fi= sj{n + uf -{Zaf. (15) 
Note that there is no singularity in this formula at = 2 because 

lim C^ix) = cos(n6) . 



III. ASYMPTOTICS 

Let us derive the Coulomb Green's function G„r(^; '^'\^) of the Schrodinger equation for 
d space dimensions. In order to do this we calculate the nonrelativistic asymptotics of the 
Coulomb Green's function of the Klein-Gordon equation valid at \E\ <^ m and Za ^ 1, 
where E = e — m. Neglecting {ZaY in /x, using the formula of summation (cf. j^), 

^0 = 5^(-l)"^C„^(a:)J2(n+.)(y) 

n=0 



23''+l/2r(z/+ 1)^^^-1/2 

and multiplying by 2m, we obtain 



7r2/2v^_i/2(w) jl + x 



irrt I I t\j 
Gnr{r,r'\E) = -j^^-^^ J ds 





sm Ks, 



X exp[2iZams + iK,{r + r') cot{Ks) — inv^ "^^^ 
K = V—^mE . 



w 



u-l/2 



(17) 



This formula is in agreement with the corresponding result of Ref.j6|. 

At high energies and small scattering angles of the particles, the characteristic angular 
momenta are large and the quasiclassical approximation is applicable. The quasiclassical 
Green's function of the Dirac equation in a Coulomb potential for d = 3 was first derived 

h [1 fi n 

in Refs. [7|, |8[. Another representation of this function was obtained in Refs. [9|, |10[. The 
quasiclassical Green's function for arbitrary spherically symmetric localized potential was 



found in Refs. 
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12| . In Ref. 13|, the quasiclassical Green's function for arbitrary localized 



potential was found with the next-to-leading quasiclassical correction taken into account. 
In Ref. jisjl . a spherical symmetry of the potential was not required. Let us consider the 
quasiclassical Green's function of the Dirac equation in a Coulomb potential for arbitrary 
space dimension d. In this case e ^ m and 1 + a; ^ 1 so that the main contribution to the 
sum over n in Eq.f lT2|) is given by ^ 1. Thus we can neglect the term {ZaY in 7, Eq.([7]), 
and perform summation over n analytically. We need to calculate two sums, 



Sa = 5^(-l)"(z/ + n + l/2)A„(x)J2(„+.+i/2)(2/), 

n=0 

00 

Sb = '^{-^TBn{x)J2(n+u+l/2){y) , 



(18) 



n=0 



where the functions An{x) and Bn{x) are defined in Eq.Q. Using the recurrent relations 
for the Bessel functions and the Gegenbauer polynomials, it is easy to show that 



[I + x)^Sb + + II2)Sb , Sb = ---^So 
ox y ox 



(19) 



so that 



S 



Try 



23'^+5/2r(i/ + l)w;'^-V2 ' 23'^+3/2r(i/ + 1)«;'^+V2 • 

Substituting these results in Eq.( [T2|) . we arrive at the final expression to the quasiclassical 



"'+'J.+,/2iw) 



(20) 



Green's function 



Ggcir,r'\e) 



1 



2!^+3/2^zy+l/2 



ds 



u 



v-1/2 



P 



sinh(ps) 



2u+l 



7 



X exp[2iZaes + ip{r + r') coth(ps) — ivrz/] Ai 

P 

M. = Ju^i/2iu) 7°£: + m — - coth(ps)(7n — 771') 



, .Ju+l/2{u) 
+ « ^ 



2 sinh [ps) 



u 



(7n + 7n') — ^0:7 pcoth(ps)[l + (an) (an')] 



p^2rr'(l + x) 

M = ■ ^ . (21) 



where p =_Y £^ — ""^^ = ii^- For d = 3 the result fl2T|) is in agreement with that obtained 
in Refs.0, [a]- The term (2'a)^ in 7, Eq.([7]), can be also neglected in the nonrelativistic 
approximation when Za ^ 1, p ^ m, and Zam/p is fixed. In this case we immediately 
obtain from Eq. (!2T1) that the nonrelativistic approximation for the Green's function of the 
Dirac equation reads 

G{r, r'\m + E) = 2:!^G„,(r, r'\E) , (22) 



where Gnrir,r'\e) is defined in Eq.([Tj 

To summarize, we have calculated in d space dimensions the Green's functions of the 
Dirac, Eq. f|T2|) . and Klein-Gordon, Eq.f fTSj) . equations in the Coulomb field. Nonrelativistic 
and quasiclassical limiting cases of these Green's functions are considered in detail. The 
results obtained can be applied for calculation of various QED amphtudes in the strong 
Coulomb field with the use of dimensional regularization. 
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